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Abstract 

Cosmic strings arising from breaking of the U{1)b-l gauge symmetry that occurs in a wide variety 
of unified models can carry zero modes of heavy Majorana neutrinos. Decaying and/or repeatedly 



> 

O \ self-interacting closed loops of these "i? — L" cosmic strings can be a non-thermal source of heavy 



right-handed Majorana neutrinos whose decay can contribute to the observed baryon asymmetry 



^fs . of the Universe (BAU) via the leptogenesis route. The B — L cosmic strings are expected in GUT 

> , 

■^ , models such as 50(10), where they can be formed at an intermediate stage of symmetry breaking 

^> ' well below the GUT scale ~ 10^^ GeV; such light strings are not excluded by the CMB anisotropy 

^: 

^^ ■ data and may well exist. We estimate the contribution oi B — L cosmic string loops to the baryon- 

>— ^ to-photon ratio of the Universe in the light of current knowledge on neutrino masses and mixings 

0^1 implied by atmospheric and solar neutrino measurements. We find that B — L cosmic string 

ij'' loops can contribute significantly to the BAU for U{1)b-l symmetry breaking scale tjb-l ~ 1-7 x 



10 GeV. At the same time, in order for the contribution of decaying B — L cosmic string loops 



k>( I not to exceed the observed baryon-to-photon ratio inferred from the recent WMAP results, the 

C^ ' lightest heavy right-handed Majorana neutrino mass Mi must satisfy the constraint Mi < 2.4 x 

10^^ {vb-l/W^^ GeV) GeV. This may have interesting implications for the associated Yukawa 
couplings in the heavy neutrino sector and consequently for the light neutrino masses generated 
through see-saw mechanism. 
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I. INTRODUCTION 



B) asymmetry of the Universe (BAU) 



A very attractive scenario of origin of the baryon 
is that it arose from a lepton (L) asymmetry jj, y, |3| • The conversion of the L-asymmetry 
to the S-asymmetry occurs via the high temperature behavior of the B + L anomaly of the 
Standard Model j^. This is an appealing route for several reasons. First, the extremely small 
neutrino masses, suggested by the atmospheric neutrino- 5], solar neutrino- 6] and Kam- 
LAND experiment |7j data, point to the possibility of Majorana masses for the neutrinos; 
such small neutrino mass can be generated, for example, through the see-saw mechanism J8| 
that involves heavy right-handed neutrinos whose interactions involve L violation in a nat- 
ural way. Second, most particle physics models incorporating the above possibility demand 
new Yukawa couplings and also possibly scalar self- couplings; these are the kind of couplings 
which, unlike gauge couplings, can naturally accommodate adequate CP violation, one of 
the necessary ingredients [9| for generating the BAU. 

Most proposals along these lines rely on out-of-equilibrium decay of the thermally gen- 
erated right-handed heavy Majorana neutrinos in the early Universe to generate the L- 
asymmetry. The simplest possibility to implement this scenario is to extend the Standard 
Model (SM) by the inclusion of a right handed neutrino, uji. A more appealing alternative 
is to consider this within the context of unified models with an embedded U{1)b-l gauge 
symmetry. For example, it can be the Left-Right symmetric model 10|, lUl where B — L is 
naturally required to be a gauge charge, or it can be a Grand Unified Theory (GUT) based 
on SO{10) gauge group. Because B — L is a gauge charge in such models, no primordial B—L 
can exist as long as the U{1)b-l gauge symmetry remains unbroken. However, spontaneous 
breaking of the U{1)b-l gauge symmetry gives heavy Majorana mass to the right-handed 
neutrinos, and a net B — L can be dynamically generated through out-of-equilibrium decay 
of these heavy right-handed Majorana neutrinos. Rapid violation of -B + L by the high 
temperature sphaleron transitions erases any B + L generated earlier. These sphaleron tran- 
sitions, however, conserve B — L. Thus, in this scenario the final BAU is related to the 
B — L produced after the U{1)b-l symmetry breaking phase transition. 

One of the interesting features of any f/(l) gauge symmetry breaking phase transition 
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12l Il3j . It has been noted 



in the early universe is the possible formation of cosmic strings 

earlier by several authors IJ, ll^ lla, ll2l, ll^ that decaying, collapsing, or repeatedly self- 



intersecting closed loops of such cosmic strings, can be a non-thermal source of massive 
particles that "constitute" the string, and that the decay of these massive particles can give 
rise to the observed BAU or at least can give significant contribution to it. Cosmic strings 
formed at a phase transition can also influence the nature of a subsequent phase transition 
that ™a. have .^pottant ^pl.cation. fot the ge.e.ation of BAU HL 

In the present context, the "5 — L" cosmic strings associated with the U{1)b-l symmetry 
breaking phase transition mentioned above are of particular interest [13 because they can 



carry zero modes [21|, |22| of the heavy right-handed neutrinos z/^. This is possible because 
the higgs field involved in the cosmic string solution arising from the spontaneous breaking 
of the U{1)b-l is the same higgs that gives heavy Majorana mass to the uji through Yukawa 
coupling. It has, therefore, been suggested [l7| that decaying closed loops of such cosmic 
strings can be an additional, non-thermal, source of the z/^j's, whose subsequent decay can 
contribute to the BAU through the leptogenesis route. 

In this paper we revisit this scenario of generating the BAU through the decay of z/^j's 
released from B — L cosmic string loops in the light of recent ideas about neutrino masses 
and mixings implied by the solar and atmospheric neutrino data. We find that, B — L 
cosmic string loops can contribute significantly to the BAU for U{1)b-l symmetry breaking 
scale TjB^L ~ 1-7 x 10^^ GeV. At the same time we show that, in order for the contribution 
of decaying B — L cosmic string loops not to exceed the observed baryon-to-photon ratio 
inferred from the recent WMAP results J22|, the lightest heavy right-handed Majorana 

1 /9 

neutrino mass Mi must satisfy the constraint Mi < 2.4 x lO"*^^ (r^^-L/lO^^ GeV) GeV. 

The rest of this paper is organized as follows. In section II we briefly discuss some exam- 
ples of symmetry breaking schemes in unified models with an embedded gauged U{X)b-l 
which potentially allow cosmic string solutions, and discuss an explicit example of a cosmic 
string solution in the context of a simple extension of the SM, namely, the gauge group 
SM ®U{1)b~l spontaneously broken to SM. The nature of the neutrino zero modes in pres- 
ence of such a cosmic string is then discussed. In section III we briefly review the evolution 
of cosmic strings with particular attention to formation of closed loops and their subse- 
quent evolution, and the production of massive particles from decaying and/or repeatedly 
self-intersecting cosmic string loops. We also discuss the observational constraints on the 
relevant cosmic string parameters. We then estimate, in section IV, the contribution of the 
B — L cosmic string loops to the BAU, and discuss the constraint on the lightest heavy right- 



handed Majorana neutrino mass Mi. Section V concludes the paper with a brief summary 
of our main results. Throughout this paper we use natural units with h = c = ks = 1- 

II. u{i)b-l cosmic strings and neutrino zero modes 

There are several realistic particle physics models where a gauged B — L symmetry exists 
and breaks at a certain scale. Since SO{10) minimally incorporates B — L gauge symmetry 
we consider the models embedded in 5*0(10). The following breaking schemes can potentially 
accommodate cosmic strings. One of the breaking schemes, motivated by supersymmetric 
50(10) J2J, |2^, involves the intermediate left-right symmetric model: 

50(10) 54 + 45 SU{?>)c®SU{2)l®SU{2)r®U{1)b-l 
126 + 126 5f/(3)c ® SU{2)l ® U{1)y O ^2 
10 + 10' SU{?,)c®U{1)q®Z2. (1) 

During the first phase of symmetry breaking, presumably at a GUT scale of ~ 10^^ GeV, 
monopoles are formed. However, during the second and third phases of symmetry breaking 
cosmic strings are formed since 7ri(^p) = 7ri(^) = Z2, where the numbers inside the 
parentheses symbolize the group structures. The monopole problem in this model can be 
solved by using a hybrid inflation ending at the left-right symmetric phase of the Universe 25| 
thus inflating away the monopoles. The formation of cosmic strings in the later phases is 
of great interest since these "light" (i.e., lighter than GUT scale) cosmic strings do not 
conflict with any cosmological observations. The Z2 strings of the low energy theory was 
investigated in an earlier work |26l |. 

Another scheme is to break supersymmetric 50(10) directly to 5f/(3)c- ® SU{2)l ® 
U{1)r ® U{1)b-l with the inclusion of extra 54 + 45. The flrst 45 acquires a vacuum 
expectation value along the direction oi B — L. However, the latter 45 acquires a vacuum 
expectation value along the direction of T^r: 

30(10) 54 + 45 + 54' + 45,' SU{3)c ® SU{2)l ®U{l)n® U{1)b-l 



126 + 126 SU{3)c®SU{2)l®U{1)y®Z2 

10 + 10' SU{3)c®U{1)q®Z2 (2) 

For the present purpose it is sufficient to consider a model based on the gauge group SM® 



U{1)b-l which is spontaneously broken to SM. Existence of cosmic strings and the related 
zero-modes in this model can be established as follows. Let the gauge field corresponding 
to the U{1)b-l symmetry be denoted by C^, and the symmetry be broken by a SM singlet 
X- Let (x) be tjb-l below the critical temperature Tb^l- In a suitable gauge a long cosmic 
string oriented along the 2;-axis can be represented (in cylindrical polar coordinates) by the 
ansatz [13| 

X = VB-Lf{r)e''^\ (3) 

where n is an integer giving the winding number of the phase of the complex higgs field Xi 
and a is the gauge coupling constant for the group U{1)b-l- In order for the solution to be 
regular at the origin we set /(O) = (7(0) = 0. Also requiring the finiteness of energy of the 
solution, we set /(r) = g{r) = 1 as r — i> 00. It turns out that both /(r) and g{r) take their 
asymptotic values everywhere outside a small region of the order of ^^i^ around the string. 
Thus away from the string (x) = ttjb-l up to a phase, and C^ is a pure gauge. The mass 
scale of the string is fixed by the energy scale of the symmetry breaking phase transition 
TjB-L at which the strings are formed. Then the mass per unit length of a cosmic string, /x, 
is of order rj^B-L ~ '^b-l- 

The lagrangian for the right-handed neutrino is 

-C^^ = iVJia'^D^UR - - [ih-JTHx^R + H.C] , (5) 

where h is the Yukawa coupling constant, a^ = {—I,a^), and i/^ = icr'^v^ defines the Dirac 
charge conjugation operation. The resulting equations of motion have been shown J2l| to 
possess \n\ normalisable zero- modes in winding number sector n. The field equations in the 
U{1) example are [27 1 

-e''[dr + '-de + ^^] d, + dt \ 

where the expressions ^ and (JH) have been substituted for x and C^, and Mr = hrjB-L- 
In the winding number sector n the normalisable zero-modes obey a^ip = tp and are of the 
form 

Mr, ^) = ( I iUir)e''' + V*{r)e'^''-'-^^') , (7) 



where U{r) and V{r) are well behaved functions at the origin and have the asymptotic 
behavior ~ exp{—Mfjr)/^/r. When nontrivial z and t dependences are included, these 
modes have solutions that depend on z + 1 and are Right movers. For n < 0, normalisable 
solutions obey a^ip = —ip, and form the zero-energy set of a Left moving spectrum. On a 
straight string these modes are massless. However on wiggly strings they are expected to 
acquire effective masses proportional to the inverse radius of the string curvature. 

III. EVOLUTION OF COSMIC STRINGS: FORMATION AND EVOLUTION OF 
CLOSED LOOPS AND PRODUCTION OF MASSIVE PARTICLES 



A. Scaling solution and closed loop formation 

The evolution of cosmic strings in the expanding Universe has been studied extensively, 
both analytically as well as numerically; for a text-book review, see the monograph (l3l |. 
Here we briefly summarize only those aspects of cosmic string evolution that are relevant for 
the subject of the present paper, namely the formation and subsequent evolution of closed 
loops of strings and production of massive particles from them. This closely follows the 
discussion in section 6.4 of Ref. |28l |. 

Immediately after their formation at a phase transition, the strings would in general be 
in a random tangled configuration. One can characterize the string configuration in terms 
of a coarse-grained length scale C,s such that the overall string energy density ps is given by 
ps = f^/^s- Initially, the strings move under a strong damping force due to friction with the 
ambient thermal plasma. In the friction dominated epoch a curved string segment of radius 
of curvature r acquires a terminal velocity oc 1/r. As a result the strings tend to straighten 
out so that the total length of the strings decreases. Thus the overall energy density in 
the form of strings decreases as the Universe expands. This in turn means that the length 
scale C,s increases. Eventually, C,s reaches the causal horizon scale ~ t. After the damping 
regime ends (when the background plasma density falls to a sufficiently low level as the 
Universe expands), the strings start to move relativistically. However, causality prevents 
the length scale ^s from exceeding the horizon size ~ t. Analytical studies supported by 
extensive numerical simulations show that the subsequent evolution of the system is such 
that the string configuration reaches a "scaling regime" in which the ratio ^ = x remains 



a constant. Numerical simulations generally find the number x to lie approximately in the 
range ~ 0.4-0.7. This is called the scaling regime because then the energy density in the 
form of strings scales as, and remains a constant fraction of, the energy density of radiation 
in the radiation dominated epoch or the energy density of matter in the matter dominated 
epoch both of which scale as t^^. 

The fundamental physical process that maintains the string network in the scaling con- 
figuration is the formation of dosed loops which are pinched off from the network whenever 
a string segment curves over into a loop, intersecting itself. In the "standard" picture |l3[, 
the closed loops so formed have average length at birth 

Lb = KTGi^t , (8) 

and they are formed at a rate (per unit volume per unit time) which, in the radiation 
dominated epoch, is given by 



dt x^ 



{TGfiy'K-H-\ (9) 



where F ~ 100 is a geometrical factor that determines the average loop length, and i^ is a 
numerical factor of order unity. 

The whole string network consisting of closed loops as well as long strands of strings 
stretched across the horizon gives rise to density fluctuations in the early Universe which 
could potentially contribute to the process of formation of structures in the Universe. More 
importantly, they would produce specific anisotropy signatures in the cosmic microwave 
background (CMB). Using a large-scale cosmic string network simulation and comparing 
the resulting prediction of CMB anisotropics with observations, a recent analysis 29| puts 
an upper limit on the fundamental cosmic string parameter /i, giving Gfi <, 0.7 x 10~^. This 
translates to an upper limit, rj ^ 1.0 x 10^^ GeV, on the symmetry-breaking energy scale of 
the cosmic string-forming phase transition. This probably rules out cosmic string formation 
at a typical GUT scale ~ 10^^ GeV. However, lighter cosmic strings arising from symmetry 
breaking at lower scales, such as the B — L cosmic strings in the case of the 5*0(10) model 
discussed in the previous section, are not ruled out. 

It should be noted here that, in the standard scenario of cosmic string evolution described 
above, the loops are formed on a length scale that is a constant fraction of the horizon length, 
as given by equation (jH}. Thus, the average size of the newly formed loops increases with 



time. At the relevant times of interest, these loops, although small in comparison to the 
horizon scale, would still be of macroscopic size in the sense that they are much larger than 
the microscopic string width scale w ~ 77^^ ~ /x^^^^. 

In contrast, results of certain Abelian Higgs (AH) model simulations of cosmic string 
evolution J30| seem to indicate that scaling configuration of the string network is maintained 
primarily by loops formed at the smallest fixed length scale in the problem, namely, on the 
scale of the width w ~ r]~^ ~ yU~^/^ of the string. These microscopic "loops" quickly decay 
into massive particles (quanta of gauge bosons, higgs bosons, heavy fermions etc.) that 
"constitute" the string. In other words, in this scenario, there is essentially no macroscopic 
loop formation at all; instead, the scaling of the string network is maintained essentially by 
massive particle radiation. In order for the scaling configuration of the string network to be 
maintained by this process, the microscopic loops must be formed at a rate 

dnb\ 1 1/2.-3 



5f = -y''t-^ ■ (10) 

The above scenario of cosmic string evolution in which massive particle radiation rather 
than gravitational radiation plays the dominant role is, however, currently a subject of 
debate [^. One of the major problems hindering a resolution of the issues involved is the 
insufficient dynamic range possible in the currently available AH model simulations and 
the consequent need for extrapolation of the simulation results to the relevant cosmological 
scales, which is not straightforward. In this paper, we shall primarily restrict ourselves to 
consideration of the "standard" macroscopic loop formation scenario described by equations 
(jHl) and above, although we shall have occasions to refer to the massive particle radiation 
scenario below (see, in particular, section III.B.2). 

B. Fate of the closed loops and massive particle production 

The behavior of the closed loops after their formation may be broadly categorized into 
following two classes: 

1. Slow death 

Any closed loop of length L in its center of momentum frame has an oscillation period 
L/2 ^]. However, a loop may be either in a self-intersecting or non-selfintersecting configu- 



ration. In general, a closed loop configuration can be represented as a superposition of waves 
consisting of various harmonics of sin's and cos's. Some explicit low harmonic number ana- 
lytical solutions of the equations of motion of closed loops representing non-selfintersecting 



loops are known in literature 
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34 l35| . and it is possible that there exists a large 
class of such non-selfintersecting solutions. Indeed, numerical simulations, while limited 
by spatial resolution, do seem to indicate that a large fraction of closed loops are born in 
non-selfintersecting configurations. 

A non-selfintersecting loop oscillates freely. As it oscillates, it loses energy by emitting 
gravitational radiation, and thereby shrinks. When the radius of the loop becomes of the 
order of its width w ~ r]^^ ~ /U^^^^j the loop decays into massive particles. Among these 
particles will be the massive gauge bosons, higgs bosons, and in the case of the B — L strings, 
massive right-handed neutrinos (z/ij) which were trapped in the string as fermion zero modes. 
We shall hereafter collectively refer to all these particles as X particles. We are, of course, 
interested here only in the urs. In addition to those directly released from the loop's final 
decay, there will also be some urs coming from the decays of the gauge and higgs bosons 
released in the final loop decay. It is difficult to calculate exactly the total number of urs so 
obtained from each loop, but we may expect that it would be a number of order unity. For 
the purpose of this paper we shall assume that each final demise of a loop yields a number 
Nfyf ~ 0(1) of heavy right handed Majorana neutrinos; we shall keep this number Nj^ as a 
free parameter in the problem. 

The rate of release of urs at any time t by the above process can be calculated as follows. 
The lifetime of a loop of length L due to energy loss through gravitational wave radiation 
is 

rGw~(rG/i)-^L. (11) 

Equations (jHl) and pi|) thus show that loops born at time t have a lifetime ~ Kt -^ H^^{t), 
where H~^{t) ~ t is the Hubble expansion time scale. It is thus a slow process. From the 
above, we see that the loops that are disappearing at any time t are the ones that were 
formed at the time {K + l)~^t. Taking into account the dilution of the number density of 
loops due to expansion of the Universe between the times of their birth and final demise, 
equation Q gives the number of loops disappearing due to this "slow death" (SD) process 



per unit time per unit volume at any time t (in the radiation dominated epoch) as 

%H = ku\ (TG,)-' (A±}l!lt-' = MK + 1)'/=^ . (12) 

at x^ K at 

where /so is the fraction of newly born loops which die through the SD process. 

The rate of release of heavy right-handed neutrinos (we shall hereafter denote them by 
A^; see section IV below) due to SD process can then be written as 

2. Quick death 

Some fraction of the loops may be born in co nfig urations with waves of high har- 
monic number. Such string loops have been shown [2a] to have a high probability of self- 



intersecting. Ref. 



gives the self-intersecting probability of a loop as 



SI 



1 - e-"-^^ , (14) 



where a = 0.4, (3 = 0.2, and N is the harmonics number. 

A self-intersecting loop would break up into two or more smaller loops. The process 
of self-intersection leaves behind "kinks" on the loops, which themselves represent high 
harmonic configurations. So, the daughter loops would also further split into smaller loops. 
If a loop does self- intersect, it must do so within its one oscillation period, since the motion 
of a loop is periodic. Under this circumstance, since smaller loops have smaller oscillation 
periods, it can be seen that a single initially large loop of length L can break up into a debris 
of tiny loops of size rj^^ (at which point they turn into the constituent massive particles) 
on a time-scale ~ L. Equation (jH} then implies that a loop born at the time t in a high 
harmonic configuration decays, due to repeated self-intersection, into massive particles on a 
time scale tqd ~ KTGfit <^ H~^(t). It is thus a "quick death" (QD) process — the loops 
die essentially instantaneously (compared to cosmological time scale) as soon as they are 
formed. Equation (0), therefore, directly gives the rate at which loops die through this quick 
death process: 

where /qd is the fraction of newly born loops that undergo QD. 



Note that, since these loops at each stage self- intersect and break up into smaller loops 
before completing one oscillation, they would lose only a negligible amount of energy in grav- 
itational radiation. Thus, almost the entire original energy of these loops would eventually 
come out in the form of massive particles. 

Assuming again, as we did in the SD case, that each segment of length ~ w ~ /i^^/^ of 
the loop yields a number N^ ~ 0{1) of heavy right-handed Majorana neutrinos, we can 
write, using equations (fT^ . Q and (jH)), the rate of release of the A^'s due to QD process as 



^ ^.„,,„J,„V.-3 



A^7v/QD:::^/i'/'t-'. (16) 

QD 



It is interesting to note here that if all loops were to die through this QD process, i.e., if 
we take /qd = 1 in equations (fTH|) and (fT?)|) . then the situation is in effect exactly equivalent 
to the microscopic loop formation scenario described by equation ()1()|1 . although the primary 
loops themselves are formed with macroscopic size given by equation (jH)). 

While the important issue of whether or not massive particle radiation plays a dominant 
role in cosmic string evolution remains to be settled, the standard model may, of course, 
still allow a small but finite fraction, /qd ^ 1, of quic kly dying loops. There already exist, 
however, rather stringent astrophysical constraints J2S!, Oil on /qd from the observed flux of 
ultrahigh cosmic rays (UHECR) above 10^^ GeV J2S| and the cosmic diffuse gamma ray back- 
ground in the energy region 10 MeV - 100 GeV measured by the EGRET experiment J39|. 
This comes about in the following way: 

The massive X particles released from the string loops would decay to SM quarks and 
leptons. The hadronization of the quarks gives rise to nucleons and pions with energy up 
to ~ Mx, the mass of the relevant X particle. The neutral pions decay to photons. These 
extremely energetic nucleons and photons, after propagating through the cosmic radiation 
background, can survive as ultrahigh energy particles. The observed flux of UHECR, there- 
fore, puts constraints on the rate of release of the massive X particles, thereby constraining 
/qd. The most stringent constraint on /qd, however, comes from the fact that the elec- 
tromagnetic component (consisting of photons and electrons/positrons) of the total energy 
injected in the Universe from the decay of the X particles initiates an electromagnetic cas- 
cade process due to interaction of the high energy electrons/positrons and photons with the 
photons of the various cosmic background radiation fields (such as the radio, the microwave 
and the infrared/optical backgrounds); see, e.g., Ref. |23| for a review. As a result, a signif- 



icant part of the total injected energy cascades down to lower energies. The measured flux 
of the cosmic gamma ray background in the 10 MeV - 100 GeV energy region J39| then puts 
the constraint [28[ 

/qd^/L < 9.6 X 10-^ (17) 

where ?7i6 = {rj/lO^^ GeV). For GUT scale cosmic strings with tjiq = 1, for example, 
the above constraint implies that /qd < 10~^, so that most loops should be in non- 
selfintersecting configurations, consistent with the standard scenario of cosmic string evo- 
lution. Note, however, that /qd is not constrained by the above considerations for cosmic 
strings formed at a scale rj ^ 3.1 x 10^^ GeV. 

In this context, it is interesting to note that there is no equivalent constraint on the corre- 
sponding parameter /sd for the slow death case from gamma ray background consideration. 
The reason is that, unlike in the QD case where the entire initial energy of a large loop 
goes into X particles, only ~ one X particle is released from a initially large loop in the 
SD case. This in turn makes the time dependence of the rate of release of massive particles 
oc t"'^ in the SD case (see equation (|T!?jl ). while it is oc t~^ in the QD case (see equation 
(I16|) ). Thus, while the SD process dominates at sufficiently early times, the QD process can 
dominate at relatively late times and can potentially contribute to the non-thermal gamma 
ray background. 

IV. CALCULATION OF BARYON ASYMMETRY 

A. Decay of heavy right-handed Majorana neutrinos and L-asymmetry 

The heavy right-handed Majorana neutrino decays to a SM lepton {£) and higgs (0) 
through the Yukawa coupling 

-Cy = fijii4>VRj + h.c. , (18) 

where fij is the Yukawa coupling matrix, and z, j = 1, 2, 3 for three fiavors. 

We shall work in a basis in which the right-handed Majorana neutrino mass matrix M 
is diagonal, M = diag(Mi, M2, M3). In this basis the right-handed Majorana neutrino is 
given by Nj = v^j ± z/^-, which satisfies N'^ = ^Nj. The standard see-saw mechanism 
then gives the corresponding light neutrino mass eigenstates ui, z/2, z/3 with masses mi, 1712, 



ms, respectively; these are mixtures of flavor eigenstates i/g, t'/i, '^r, and are also Majorana 
neutrinos, i.e., z/j = z/f. 

The decays of the heavy right-handed Majorana neutrinos can create a non-zero L — 
asymmetry (which is ultimately converted to 5-asymmetry) only if their decay violates CP. 
The CP-asymmetry parameter in the decay of Nj is defined as 

_ r(iv, -. £0) - r(iv, -. t^^) 

Assuming a mass hierarchy in the heavy neutrino sector. Mi < M2 < M3, it is reasonable 
to expect that the final lepton asymmetry is produced mainly by the decay of the lightest 
right handed neutrino A^^i. Any asymmetry produced by the the decay of N2 and A^3 will be 
washed out by the lepton number violating interactions mediated by the Ni. As the Universe 
expands, the temperature of the thermal plasma falls. Below a temperature Tp ~ Mi, all 
L- violating scatterings mediated by A^i freeze out, thus providing the out-of-equilibrium 
situation J2| necessary for the survival of any net L-asymmetry generated by the decay 
of the A^i's. The final L-asymmetry is, therefore, given essentially by the CP asymmetry 
parameter ei. 

An accurate calculation of the net L-asymmetry can only be done by numerically solving 
the full Boltzmann equation that includes all lepton number violating interactions involving 
all the Nj^s present at any time, including the Nj^s of non-thermal origin such as the ones 
produced from the decaying cosmic string loops, as well as those of thermal origin. This 
is beyond the scope of the present paper; here we shall simply assume that below the 
temperature Tp = Mi, all interactions except the decay of the A^i are unimportant, so that 
each A^i released from cosmic strings additively produces a net L-asymmetry ei when it 
decays. 

To fix the value of ei, we note that there is an upper bound [4(1] on ei, which is related 
to the properties of the light neutrino masses. In a standard hierarchical neutrino mass 
scenario with m.^ ^ m,2 > m,i , this upper limit is given by J40|, |41 1 



, , 3 MiTTlq , ^ 

|ei|<- V^, 20 

where v ~ 174 GeV is the electroweak symmetry breaking scale. Furthermore, the above 



upper limit is in fact saturated 
we shall assume to be the case. 



4l| in most of the reasonable neutrino mass models, which 



The atmospheric neutrino data J5[ indicate u^ ^^ z/^ oscillations with nearly maximal 
mixing (6'atm — 45°) and a mass-squared-difference Aml^^ = |m| — m^] ~ 2.6 x 10~^eV^. 
The solar neutrino [6[ and KaniLAND [7| data, on the other hand, can be explained by 
^e ^^ ^fi oscillations with large mixing angle (LMA) (6*801 — 32°) and Amg^j = \m2 — m1\ ^ 
7.13 X 10~^eV^. Assuming, again, the standard light neutrino mass hierarchy, the above 
numbers give rris ^ (Aml^^^) ~ 0.05 eV. In our calculations below, we shall use 

e, ^ 9.86 X 10- ( -^] (i^a^ 1 . (21) 



^1013 GeVy I 0.05 eV 
The L-asymmetry is partially converted to a 5-asymmetry by the rapid nonperturbative 
sphaleron transitions which violate B + L but preserve B — L. Assuming that sphaleron 
transitions are ineffective at temperatures below the electroweak transition temperature 
(Tew); the 5-asymmetry is related to L-asymmetry by the relation |42i | 

B=p(B-L) = ^^Lc^-0.55L, (22) 

p — I 



where we have taken p = 28/79 appropriate for the particle content in SM J42]. If sphaleron 
transitions continue to be effective below Tew? then the above relation between B and L is 
slightly modified; we can however ignore this at the level of accuracy aimed at in the present 
paper. 

The net baryon asymmetry of the Universe is defined as 



where 



Ys = ^^^^^^, (23) 

s 



s ~ {27ryA5)gX' ^ 43.86 (^) T' (24) 



is the entropy density, g^ being the number of relativistic degrees of freedom contributing 
to the entropy at the temperature T. At temperatures in the early Universe relevant for the 
process of baryon asymmetry generation, g^, ^ 100 in SM. 

Observationally, the BAU is often expressed in terms of the baryon-to-photon ratio 77 = 
{ub — n§)/nj, whose present-day-value t]q is related to that of Yb through the relation 

no ~ 7.0^^,0 • (25) 

The observed value of rjo inferred from the Wilkinson Microwave Anisotropy Probe (WMAP) 
data is 23[ 

^WMAP ^ ^g_^+o.3^ X 10-10. (26) 



We now proceed to estimate the contribution to the BAU from the two cosmic string 
loop processes discussed in the previous section. 

B. Slow death case 

The contribution of the SD process to rjo can be written as 

ryo^° ^ 7.0 X 0.55 eiTif^) dt , (27) 

where tp is the cosmic time corresponding to the temperature Tp ^ Mi and to is the present 
age of the Universe. Using equations (fT^ . P^ and the standard time-temperature relation 
in the early Universe, 

tc.0.39.-/^^, (28) 

where Mp\ ~ 1.22 x 10^^ GeV is the Planck mass, we see that the dominant contribution to 
the integral in equation ()27|) comes from the time tp -C to, i.e., from the epoch of temperature 
Tp ^ Ml, giving 

= 2.0 X 10-' N^hf {j^^^Y , (29) 

where we have defined the Yukawa coupling hi = Mi/tje-l, used equation (^1]) for ei with 
(Am^^-jjj) = 0.05 eV, and also taken x = 0.5, T = 100, K = 1 and /sd = 1 in equation 

m- 

The Yukawa couplings are generally thought to be less than unity. With hi < 1, we see 
from (PU|) and (PH|) that the cosmic string loop slow death process can produce the observed 
BAU only for B — L phase transition scale 

Vf-L ^ 5.5 X lO^^iV"'/' GeV . (30) 

Assuming Nj^ <, 10, say, we see that cosmic string loop SD process can contribute to BAU 
if rjB~L ~ 1-7 X 10^^ GeV; lower values of tjb-l are relevant only if we allow hi > 1. 

At the same time, for a given tjb-l satisfying (|30|). in order that the contribution (|29|) not 
exceed the highest allowed observed value of tjq given by equation (|^ , the Yukawa coupling 
hi must satisfy the constraint 

1^1 ~^-^^Uoi3GeVy ^ ' ^^^ 



which, in terms of the hghtest heavy right handed Majorana neutrino mass Mi, reads 

Mf ° < 2.4 X 10-iV^^/^ (l^f^) ''' GeV . (32) 

Note the rather weak dependence of the above constraints on N^. Also, the 4th power 
dependence on Mi of equation ()29p and the rather narrow range of the observed value of 
rjQ given by equation (J26|l together imply that, in order for the SD process to explain the 
observed BAU, Mi (and equivalently hi) cannot be much smaller than their respective values 
saturating the above constraints. 

C. Quick death case 

Replacing (^^)gj-) in equation (J2Z|) by {^^)qy) given by equation (fT?)|l . and following 
the same steps as in the SD case above, we get the contribution of the QD process to tjq as 

V^ ^ 5.17 X lO-^^N^f^,, {j,^Y {j^§^) 

= 5.17 X IO-i^at^/qd/i? {j^i^f . (33) 

From (p^ and (f^H|) we see that, considering the most optimistic situation with /qd = 1, the 
QD process is relevant for BAU only for 

77^°^ > 1.1x10^4 AT^'/^^GeV; (34) 

lower values of rjB-L are relevant only if we allow hi > 1. On the other hand, the constraint 
(fTTj) allows /qd = 1 only if tjb-l < 3.1 x 10^^ GeV. This can be reconciled with the above 
constraint (jMj) only for N^ > 45 or so. Such a large value of N^ seems unlikely. 
In general, using the constraint (fTTj) on /qd in (jHHj) we get 

= 5.0 X IQ-^^N^hl (it^g^) . (35) 

Comparing again with the observed value of r/o, we now see that, for hi < 1, the QD process 
can be relevant for BAU only for 

rff_^ > 1.2 X lO'^iV^^ GeV . (36) 



For values of tjb-l satisfying the above constraint (j^ . the QD process can produce the 
observed value of BAU for 

which in terms of Mi now reads 

Mr ^ 3.6 X lO^^iV^^/^ (l^f^) ''' GeV . (38) 

From the above discussions we see that, as far as their contributions to the BAU is con- 
cerned, the QD process becomes important only at relatively higher values of the symmetry 
breaking scale tje-l compared to the SD process. 

V. SUMMARY AND CONCLUSIONS 

A wide class of unified theories with an embedded U{1)b~^l gauge symmetry allows for- 
mation of "-B — L" cosmic strings at the U{1)b-l symmetry-breaking phase transition at a 
symmetry-breaking scale tjb-l well below the GUT scale (~ 10^^ GeV). Such "light" cosmic 
strings are not currently excluded by CMB anisotropy data. The B — L cosmic strings can 
carry zero modes of heavy right handed Majorana neutrinos (A^), and the latter can be 
released from closed loops of these cosmic strings when the loops eventually disappear. The 
decay of the A^'s can give rise to a L-asymmetry which is partially converted to 5-asymmetry 
via nonperturbative sphaleron transitions. 

In this paper we have studied the contribution to the baryon asymmetry of the Universe 
due to decay of heavy right handed Majorana neutrinos released from closed loops of i? — L 
cosmic strings in the light of current ideas on light neutrino masses and mixings implied 
by atmospheric and solar neutrino measurements. We have estimated the contribution to 
BAU from cosmic string loops which disappear through the process of (a) slow shrinkage 
due to energy loss through gravitational radiation — which we call slow death (SD), and (b) 
repeated self-intersections — which we call quick death (QD). We find that for reasonable 
values of the relevant parameters, the SD process dominates over the QD process as far 
as their contribution to BAU is concerned. We find that B — L cosmic string loop SD 
process can contribute significantly to, and can in principle produce, the observed BAU for 
U{1)b-l symmetry breaking scale tjb-l ~ 1-7 x 10^^ GeV. The QD process, on the other 
hand, becomes relevant for BAU only for relatively higher values of tjb-l ~ 10^^ GeV. 



We have also found that, in order for the contribution of decaying B — L cosmic string 
loop SD process not to exceed the observed baryon-to-photon ratio inferred from the recent 
WMAP results, the lightest heavy right-handed Majorana neutrino mass Mi must satisfy the 
constraint Mi < 2.4 x 10^^ ('7_b-l/10^'^ GeV) GeV. This result may have interesting impli- 
cations for the associated Yukawa couplings in the heavy neutrino sector and consequently 
for the light neutrino masses generated through see-saw mechanism. 

We conclude that processes involving closed loops of cosmic strings formed at a U{\)b-l 
symmetry breaking phase transition may make significant contribution to the observed BAU 
and should be included in considerations of baryon generation processes in general. A full 
Boltzmann equation calculation of the baryogenesis process including the heavy right handed 
neutrinos of cosmic string origin as well of thermal origin is in progress and will be reported 
elsewhere. 
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